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Introduction 

The  flows  of  a  weakly  non-equilibrium  compressible  gas  are  satisfactorily  described  by 
the  Navier  -  Stokes  equations.  The  stress  tensor  in  these  equations  contains  two  dissipative 
coefficients  that  are  dynamic  viscosity  and  bulk  viscosity  //2  .  The  last  one  takes  into 
account  relaxation  of  the  excited  internal  degrees  of  freedom  of  gas  molecules. 

In  paper  [1]  some  experimental  results  on  the  bulk  viscosity  influence  on  the  value  of  the 
critical  Reynolds  numbers  of  the  laminar-turbulent  transition  (LTT)  are  presented.  From  the 
experimental  data  obtained  it  follows  that,  with  increasing  bulk  viscosity,  the  value  Rec  of 
LTT  also  grows.  In  this  case,  registered  change  reached  ten  percent.  However,  the  matter  of  the 
reliability  of  the  obtained  data  have  been  debatable.  The  calculations,  which  were  derived 
within  the  framework  of  linear  stability  theory  in  papers  [2,  3],  have  shown  that  bulk  viscosity 
influence  on  LTT  is  rather  little.  In  [2]  the  calculations  were  carried  out  for  air  at  Mach  number 
M0  =  4;  5  and  relation  of  bulk  viscosity  to  dynamic  viscosity  a  =  //2  /  >]\  didn’t  exceed  unit.  In 
paper  [3]  the  estimations  were  obtained  up  to  a  <  30  . 

At  the  same  time  the  data  of  the  experiments  [1]  correspond  immediately  to  the  nonlinear 
stage  of  transition  in  the  turbulent  regime  at  a  «7  .  It’s  known  from  the  modern  scenarios  of 
the  transition  and  generation  of  turbulence  are  that  these  processes  are  essential  nonlinear.  They 
realize  by  means  of  origin,  evolution  and  decay  of  certain  vortex  structures  [4].  In  boundary 
layers  and  pipes  they  are  known  X  -structures,  in  two-dimension  free  shear  flows  and  jets  they 
are  two-dimensional  vortices  extended  on  the  transversal  coordinate.  Under  this  interpretation 
the  generation  of  turbulence  can  be  considered  as  the  process  of  LTT  which  is  stochastically 
repeated  in  space  and  time.  This  fact  allows  one  to  assume  that  bulk  viscosity  influence  on  the 
development  of  perturbations  at  the  nonlinear  stage  can  be  estimated  by  simulating  the 
interaction  of  the  solitary  vortex  structure  with  the  basic  (mean)  flow.  In  this  paper  the  simple 
model  of  the  evolution  of  the  spanwise  vortex  structure  in  a  shear  flow  of  non-equilibrium 
molecular  gas  is  considered  for  such  a  purpose. 

Flow  model  and  statement  of  problem 

The  modeling  problem  is  posed  as  follows.  The  gas  shear  flow  is  considered  in  a  two- 
dimension  square  area  with  the  side  / .  The  unperturbed  velocity  field  is 

U  =  2U0yi/l,  y  e[-//2;//2]. 

Here  i  is  a  unit  vector  along  axis  x  ,  U0  is  an  absolute  value  of  a  velocity  vector  at  the  upper 
and  lower  bounds  of  a  cell.  The  finite  amplitude  of  the  perturbation,  which  is  a  circle  cross 
vortex  with  radius  RQ  with  uniform  vorticity  Q(|  centered  in  the  origin  of  coordinates  is 

imposed  on  the  basic  flow  at  the  moment  tQ.  Then  the  evolution  of  the  structure  in  a 
compressible  viscous  flow  of  thermal  exited  gas  is  considered. 
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Under  an  assumption  of  low  excitation,  the  evolution  of  the  structure  in  a  modeling  cell  is 
described  by  the  system  of  the  Navier  -  Stokes  equations  of  compressible  gas.  It  is  assumed  that 
the  transport  coefficients  are  constant.  In  dimensionless  variables 

t  =7  /  t0,x  =x  /2R0,y  =y  /2R0,ux  =ux  IU0,uy  =  uy  /U0, 

P  =  P!P0T  =  T  /T0,p  =  ~p / pQU q,t0  =2 R0/U0 

the  system  of  the  equations  is  written  as  follows 
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Here  MQ=UQ/c0  is  the  Mach  number,  Re  =  2 p0UQRQ/t]1  is  the  Reynolds  number, 
Pr  =  t)lcp  /  Aq  is  the  Prandtl  number,  cQ  is  the  adiabatic  sonic  speed,  c  and  AQ  are  the 

specific  heat  at  constant  pressure  and  the  conductivity  respectively,  and  y  is  the  ratio  of 

specific  heats.  The  initial  -  boundary  conditions  for  system  (1)  were  chosen  in  the  following 

1 

form.  At  the  initial  moment,  the  components  of  the  vector  velocity  u ,  temperature  T  and 
density  p  are  correspondingly  equal 
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Here  the  parameter  /3  =  Q0R0/2Uq  characterizes  the  intensity  of  the  imposed  vortex 

perturbation,  and  the  relation  A  =  l  I2RQ  can  be  considered  as  an  intermittence  coefficient.  The 

periodicity  conditions  were  stated  on  the  bounds  of  a  modeling  cell  for  normal  velocity 
components,  density  and  temperature.  The  conditions  of  antisymmetry  were  used  for  tangential 
velocity  components. 

The  equation  for  integral  production  of  turbulent  energy  in  a  modeling  cell  was  used  for 
estimation  of  bulk  viscosity  influence.  It  was  deduced  on  the  analogy  with  a  similar  equation  for 
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unbounded  flows  [5].  The  velocity  fluctuations  were  defined  as  the  difference  between 
instantaneous  velocity  values  and  velocity  field  of  the  carried  gas  flow 


-  2  y  , 

Ux  =  ux  -  —  ,uy  =  uy 

and  the  pulsations  of  density  and  pressure  were  defined  by  the  relations 

i  '  1 
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For  fluctuating  variables  from  system  (1),  the  continuity  and  momentum  equations  were 
deduced,  which  are  used  to  obtain  the  required  equation  for  turbulent  energy  production  by 
integration  over  the  cell.  It  has  the  form 
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The  integrals  on  the  right-hand  side  of  the  Eq.  (2)  describe  the  interchange  of  energy  between 
the  structure  and  mean  flow,  the  work  under  fluctuating  compression  and  expansion  of  a  gas, 
and  the  dissipation  of  perturbation  energy,  respectively.  Let  us  notice  that  the  integrals  along  the 
cell  bounds  arising  during  the  derivation  of  Eq.  (2)  with  used  boundary  conditions  are  precisely 
nullified.  The  values  of  the  integrals  J {  and  J ,  can  be  either  positive  or  negative.  The  integrals 


J3  and  J 4  are  positively  defined.  From  Eq.  (2)  it  follows  that  bulk  viscosity  brings  in 

fluctuation  energy  production  a  term  of  the  same  sign  as  one  caused  by  dynamic  viscosity,  and 
both  terms  reduce  the  energy  of  the  perturbation.  A  similar  equation  was  used  for  the  energy 
estimation  of  the  critical  Reynolds  number  Re  f  of  LTT  [6]  which  was  obtained  by 
minimization  of  the  functional 
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at  dE  /  dt  =  0  . 

Therefore,  it  is  possible  to  conclude  that  Re  grows  with  an  increase  in  the  parameter  a 


(or  bulk  viscosity  77, ). 


The  method  of  the  solution  and  the  results  of  the  calculations 

For  the  numerical  solving  of  the  modeling  problem  an  explicit  -  implicit  finite  - 
difference  scheme  with  splitting  on  directions  and  physical  processes  was  used  [7].  The 
calculations  were  carried  out  for  the  Mach  number  M0  =  0.6  ,  the  Prantdl  number  Pr  =1.0  and 
the  values  of  the  parameter  a  =  0;  0.5;  1;  1.5;  2  .  The  characteristics  of  the  considered  structure 
were  parametrized  according  to  the  data  calculated  for  the  mixing  layer  [8].  The  Reynolds 
number,  the  intermittence  coefficient,  and  the  relative  intensity  of  the  vortex  structure  were 
chosen  as  Re  =  60, X  =  3,j3  =  1.4  respectively.  The  mesh  domain  contained  31x31  =  961 

nodes  with  a  step  h  =  0.1  on  both  coordinates,  the  time  step  was  r  =  0.01 .  The  evolution  was 
traced  up  to  600  time  steps. 
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Figure  1. 


To  test  the  numerical 
scheme,  diffusion  of  a  circular 
isochoric  Rankine  vortex  cen¬ 
tered  in  the  origin  of  coordi¬ 
nates  in  the  absence  of  a  shear 
flow  were  computed.  Then  the 
results  was  compared  with  the 
corresponding  analytical  solu¬ 
tion.  It  was  shown  that  the 
maximal  relative  error  did  not 
exceed  0.6%  . 

The  average  Reynolds 
stresses  on  a  cell  and  on  a  time 
interval  8  were  calculated  for 
the  estimation  of  bulk  viscos¬ 
ity  influence  on  turbulence 
characteristics 
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The  dependences  of  the  average  Reynolds  stresses  on  the  parameter  a  (or  bulk  viscosity  ijn ) 
are  presented  in  Fig.  1  for  Re  =  60,  A  =  3,  /?  =  1.4  and  8  =  6.  They  show  that  average  Reynolds 
stresses  a  decrease  if  the  value  of  bulk  viscosity  rj2  grows.  The  relative  deviation  in 
calculated  interval  a  varies  in  the  limits  0.4%  -2.7%  compared  with  a  =  0  . 

As  known,  the  organized  vortices  bring  about  40%  of  contribution  to  the  total  Reynolds 
stresses  in  the  mixing  layer  [9].  At  the  same  time,  viscous  dissipation  influence  on  large  -  scale 
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Figure  2 
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vortices  is  essentially  less  than  the  small-scale  component  of  turbulence.  This  fact  allow  us  to 
confirm  that  the  decrease  in  the  total  Reynolds  stresses  can  be  estimated  at  least  in  the  same 
limits  for  0  <  a  <  2  . 

In  Fig.  2  the  time  dependences  of  turbulent  energy  production  for  various  values  of  the 
parameter  a  derived  from  Eq.  (2)  are  plotted.  Analyzing  the  curves,  we  see  that  the  production 
of  turbulent  energy  is  always  negative  and  decreases  in  this  modeling  problem.  Moreover,  the 
higher  the  value  of  bulk  viscosity  the  more  intense  the  suppression  of  fluctuations.  The  absence 
of  the  mechanism  of  the  positive  production  of  perturbation  energy  such  as  the  stretching  of 

quasi-streamwise  vortices  in  the  mixing  layer  [9]  is  one  disadvantage  of  the  model  and  it  does 
not  allow  direct  estimation  of  bulk  viscosity  influence  on  Re  of  LTT. 

Conclusion 

The  results  of  numerical  modeling  of  vortex  structure  dynamics  in  a  shear  flow  allow  us 
to  conclude  that  bulk  viscosity  stabilizes  LTT.  It  is  shown  that  the  Reynolds  stresses  averaged 
over  space  and  time  decrease  if  bulk  viscosity  increases  with  ranging  from  r/2  =  0  up  to 

77,  =2rjl .  The  obtained  change  in  the  Reynolds  stresses  varies  from  0.4%  up  to  2.7%.  The 

production  of  turbulent  energy  in  the  considered  modeling  problem  is  always  negative,  and  the 
higher  the  value  a ,  the  greater  the  decrease  in  turbulent  fluctuations.  The  change  is  essential 
because  by  the  order  of  the  value  it  is  comparative  with  the  effect  of  practically  employed 
mechanical  devices  for  reducing  turbulent  drag,  for  example,  by  means  ofriblets  [4].  Therefore, 
the  hypothetical  opportunity  of  processing  turbulent  drag  in  compressible  flows  of  molecular 
gases  by  changing  bulk  viscosity  has  appeared.  In  connection  with  that,  it  is  necessary  to 
continue  researches  on  the  basis  of  more  improved  models. 

The  work  is  partly  financial  supported  by  the  Russian  Foundation  of  Basic  Research 
(Grant  No.  01-01-00827). 
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